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Abstract. Resonant scattering of plane waves by a periodic slab under conditions close 
to those that support a guided mode is accompanied by sharp transmission anomalies. 
For two-dimensional structures, we establish sufficient conditions, involving structural 
symmetry, under which these anomalies attain total transmission and total reflection at 
frequencies separated by an arbitrarily small amount. The loci of total reflection and 
total transmission are real-analytic curves in frequency-wavenumber space that intersect 
quadratically at a single point corresponding to the guided mode. A single anomaly or 
multiple anomalies can be excited by the interaction with a single guided mode. 
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1 Introduction 

A dielectric slab with periodically varying structure can act both as a guide of electromagnetic 
waves as well as an obstacle that diffracts plane waves striking it transversely. Ideal guided modes 
exponentially confined to a slab are inherently nonrobust objects, tending to become leaky, or to 
radiate energy, under perturbations of the structure or wavevector. The leaky modes are manifest 
as high-amplitude resonant fields in the structure that are excited during the scattering of a plane 
wave. One may think of this phenomenon loosely as the resonant interaction between guided 
modes of the slab and plane waves originating from exterior sources. This interaction generates 
sharp anomalies in the transmittance, that is, in the fraction of energy of a plane wave transmitted 
across the slab as a function of frequency and wavevector. 

In this work, we analyze a specific feature of transmission resonances for two-dimensional lossless 
periodic structures (Fig. [T]) that results from perturbation of the wavenumber from that of a true 
(exponentially confined to the structure) guided mode. Graphs of transmittance vs. frequency u 
and wavenumber k parallel to the slab typically exhibit a sharp peak and dip near the parameters 
{kq, uq) of the guided mode. Often in computations these extreme points appear to reach 100% and 
0% transmission, which means that, between two closely spaced frequencies, the slab transitions 
from being completely transparent to being completely opaque (Fig. [2]). The main result, presented 



in Theorem 12, is a proof that, if the slab is symmetric with respect to a line parallel to it (the 
j;-axis in Fig. [T]), then these extreme values are in fact attained. Subject to technical conditions 
discussed later on, it can be paraphrased like this: 

Theorem. Consider a two-dimensional lossless periodic slab that is symmetric about a line parallel 
to it. If the slab admits a guided mode at an isolated wavenumber- frequency pair {kq^loq), then total 
transmission and total reflection are achieved at nearby frequencies whose difference tends to zero 
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as K tends to kq. The loci in real {K,,u;)-space of total transmission and total reflection are real 
analytic curves that intersect tangentially at [kq^ujq). 

In this Theorem, it is important that the slab admits a guided mode at an isolated pair (kq, ooq). 
The frequency ojq is above cutoff, meaning that it hes above the hght cone in the first Brillouin 
zone of wavenumbers k and therefore in the (k, w) regime of scattering states (Fig. [s]). Perturbing 
K from destroys the guided mode and causes resonant scattering of plane waves at frequencies 
near ujq. In the literature, one encounters these nonrobust guided modes at kq = 0, that is, they 
are standing waves. Although we are not aware of truly traveling guided modes {k / 0) above 
cutoff in periodic photonic slabs, we believe that they should exist in anisotropic structures. 

Resonant transmission anomalies are well known in a wide variety of applications in electro- 
magnetics and other instances of wave propagation, and a veritable plenitude of models and tech- 
niques has been developed for describing and predicting them [3 |2l |6l |8j . The causes of anomalies 
are manifold and include Fabry-Perot resonance and Wood anomalies near cutoff frequencies of 
Rayleigh-Bloch diffracted waves. The present study addresses the specific resonant phenomenon 
associated with the interaction of plane waves with a guided mode of a slab structure. 

In our point of view, one begins with the equations of electromagnetism (or acoustics, etc.), 
admitting no phenomenological assumptions that cannot be proved from them, and seeks to provide 
rigorous theorems on the phenomenon of resonant scattering. A rigorous asymptotic formula for 
transmission anomalies in the case of perturbation of the angle of incidence (or Bloch wavevector) 
has been obtained through singular complex analytic perturbation of the scattering problem about 
the frequency and wavenumber of the guided mode for two-dimensional structures |16l [T2l I13j . 



and the analysis in this paper is based on that work. The essential new result is Theorem 12 
on total refiection and transmission. Previous analyticity results were based on boundary-integral 
formulations of the scattering problem, which are suitable for piecewise constant scatterers. Here, 
we deal with general positive, coercive, bounded coefficients and thus give self-contained proofs of 
analyticity of a scattering operator based on a variational formulation of the scattering problem. 

There are interesting open questions concerning the detailed nature of transmission resonances. 
In passing from two-dimensional slabs (with one direction of periodicity) to three-dimensional slabs 
(with two directions of periodicity), both the additional dimension of the wavevector parallel to the 
slab as well as various modes of polarization of the incident field impart considerable complexity to 
the guided-mode structure of the slab and its interaction with plane waves. The role of structural 
perturbations is a mechanism for initiating coupling between guided modes and radiation [3] jH 
§4.4] that also deserves a rigorous mathematical treatment. A practical understanding of the 
correspondence between structural parameters and salient features of transmission anomalies, such 
as central frequency and width, would be valuable in applications. 

The main theorem is proved in Section [3] in the simplest case in which the transmittance graph 
exhibits a single sharp peak and dip. The proof rests on the complex analyticity with respect to 
frequency and wavenumber inherent in the problem of scattering of harmonic fields by a periodic 
slab. The framework for our analysis is the variational (or weak-PDE) formulation of the scattering 
problem, which is reviewed in Section [2j Section |4] deals briefly with non-generic cases in which 
degenerate or multiple anomalies emanate from a single guided-mode frequency. A number of 
graphs of transmittance in the generic and non-generic cases are shown in Section [5] 
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Figure 1: An example of a two-dimensional periodic slab. One period truncated to the rectangle [— vr, tt] x L, L] 
is denoted by Q. 
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Figure 2: Numerical computation of the transmittance, or the percentage of energy transmitted across a penetrable 
waveguide of period 2-k as a function of the frequency of the incident plane wave. Here, the wavenumber in the 
x-direction (Fig. [T| is k = 0.02 and one period consists of a single circle of radius tt/2 with e = 10 and an ambient 
medium with e = 1; ^ — 1 throughout. The structure supports guided modes at = (0,0.5039...) and (K,ti;) = 

(0, 0.7452...), both contained within the region V of one propagating diffractive order (Fig[3|. Theorem 1 1 2| guarantees 
that the transmittance attains minimal and maximal values of 0% and 100% at each of the sharp anomalies near the 
guided-mode frequencies. This example was computed using the boundary-integral method described in [4]. 



2 Background: Scattering and guided modes 

Readers familiar with variational formulations of scattering problems can easily skim this section 
and proceed to Section [3j which contains the main result. 

A two-dimensional periodic dielectric slab (Fig. [T]) is characterized by two coefficients e{x,z) 
and fx{x, z), {x, z) G M^, that are periodic in the x-direction and constant outside of a strip parallel 
to the X axis. We take these coefficients to be bounded from below and above by positive numbers: 

e{x + 27rn, z) = e(x, z), -|- 27m, z) = fi{x, z), for all n G Z, 

e{x, z) = eo, ^(x, z) = fio, if \z\ > L, (1) 

< e_ < e(x, z) < e+, < ^_ < /u(x, z) < 

Physically, the structure is three-dimensional but invariant in the y-direction. The Maxwell 
system for y-independent electromagnetic fields in such a structure decouples into two polarizations 
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and simplifies to the scalar wave equation for the out-of-plane components of the E field and the 
H field independently. We will consider harmonic fields, whose circular frequency oj will always be 
taken to be positive. Given a frequency w, plane waves and guided modes are characterized by their 
propagation constant n in the direction parallel to the slab. The y component Ey of a harmonic 
i?-polarized field with propagation constant k is of the pseudoperiodic form 

Ey{x,z,t) = n(x,z)e*('^^-'^*), 

u{x + 27rn, z) = u{x, z) for all n G Z, 

in which the periodic factor u satisfies the equation 

(V + in) • (V + in)u{x, z) + e u?u{x, z) = 0, (2) 

with K = (k,0). The number n can be restricted to lie in the Brillouin zone [—Via, V2)- As long as 
the quantities 

fill = eo^ow^ - (m + Kf (3) 

are nonzero for all integers m, the general solution of this equation admits a Fourier expansion on 
each side of the slab, 

oo 

u{x, z)= {At^e'"^""' + S^e-^''™")e*'"^' for ±z> L, (4) 

m=--oo 

For real k and a; > 0, the square root is chosen with a branch cut on the negative imaginary axis, 
and the sign is taken such that rjm = \i]m\ Vm > ^ ^^'^ Vm = i\r]m\ < ^■ 

2.1 Scattering and guided modes in periodic slabs 

The theory of guided modes underlies the analysis of resonance in Section [s} We present the 
pertinent elements of this theory here and refer the reader to more in-depth discussions in [13] . 

In the problem of scattering of plane waves for real {k,uj), one takes in Q A~ = = 
for the infinite set of m such that r]"^ < to exclude fields that grow exponentially as \z\ — )• oo. 
The exponentially decaying Fourier harmonics are known as the evanescent diffraction orders. The 
finitely many propagating diffractive orders (r/^ > 0) express the sum of the incident plane waves 
and the scattered field far from the slab. In view of the factor e"*"^* and the convention cj > 0, we 
see that A~ and are the coefficients of inward traveling plane waves and and are the 
coefficients of outward traveling waves. The linear orders ijm = correspond to "grazing incidence" , 
and will not play a role in the present study. 

If for a given pair {k,u;), r]m 7^ for all m, the numbers ijm can be continued as analytic 
functions in a complex neighborhood of {k,uj) G C^. The following outgoing condition is central 
to the mathematical formulation of the scattering problem and the definition of generalized guided 
modes. 

Condition 1 (Outgoing Condition). A pseudo-periodic function u{x,z) = u[x,z)e^'^^ is said to 
satisfy the outgoing condition for the complex pair {k,uj), with Re(w) > if there exist a real 
number L and complex coefficients {c^} such that 

u{x, z) = Y, c^e^^'^™"e^™^ for ± z > L. 
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Figure 3: The diamond T> of one propagating diflractive order (m = 0) within the first BriUouin zone is shown in 
grey. The point {hio,ujo) represents an isolated real pair on the complex dispersion relation for generalized guided 
modes. If the periodic slab is symmetric about a line parallel to it, as shown in Fig.[lj the main theorem (Theorem 12 
Sec.[3| guarantees that, under certain generic conditions, the zero sets of transmission, T{k,uj) = (or b{K,Lu) = 0), 
and reflection, R{hi,uj) = (or a{n,ui) = 0), intersect quadratically at the pair (Ko,ajo), where T and R are the 
complex transmission and reflection coefficients. Since |rp + \R\^ = 1, one sees that |r|, being continuous in a 
punctured neighborhood of {k.o,ujo), achieves all values between and 1 in each neighborhood of (kqi'^o)- Each zero 
set either intersects the boundary of I?, as depicted here, or attains an infinite slope. 



We will be concerned with the case of exactly one propagating harmonic m = 0. This regime 
corresponds to real pairs {k,uj) that lie in the diamond 

P = Lj) G : |k| < 72 and \k\ < oj^equq < 1 — |k|} 

shown in Fig. [sj The numbers rj^a are analytic functions of in a complex neighborhood V of 

V; thus, M2 D 2? c P' C C^. 

The problem of scattering of plane waves by a periodic slab is the following. 

Problem 2 (Scattering problem). Find a function u{x,z) such that 

u{x,z) = e^'^^u{x, z), u is 2TT-periodic in x, 
< {V + iK)-n~'^{V + iK)u{x,z) + eio^u{x,z) = 0, for {x, z) e R'^ , (5) 
u{x,z) = u™'^{x,z) + u^'^{x, z), u^^ satisfies the outgoing condition, 

in which z) = Aqc^'"'^ + Boe-'^^o^ 

The weak formulation of the scattering problem is posed in the truncated period 

Q, = {{x, z) : —IT < x < TT, —L < z < L}, 

T = \^T±, T± = {(x, z) : —TT < X < TT, z = ±L}. 

and makes use of the Dirichlet-to-Neumann map N = N{k, uj) on the right and left boundaries 
F-i- that characterizes outgoing fields. This is a bounded linear operator from Hi(r) to H'i{T) 
defined as follows. For any / G ff2(r), let fm = ifm^fm) be the Fourier coefficients of /, that is, 
/(±L,x) = E^/±e^-^ Then 

N:H"HT)^H-'2{T), (Ff)^ = -ir?^/^. (6) 
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This operator has the property that 

dnU + Nu = on r <;=^ u is outgoing, 
where "on F" refers to the traces of u and its normal derivative on F. In the periodic Sobolev space 

H^^^{Q) = {n G H\n) : n(7r, z) = u{-7r, z) for all z E (-L, L)}, 
in which evaluation on the boundaries of 0, is in the sense of the trace, define the forms 

h{u,v) = hfi^^{u,v) = / + iK)u- {V -iK)v + j {Nu)v, (7) 

Jo, Jr 

b{u, v) = e uv, 
Jn 

P{V) = PnAy) = ^^o' I idnv""" + Nu''''')v. 



Remark on notation. The forms h and p depend on k and oj through = N{k,,u;). In the sequel, 
the dependence on k and lv of certain objects such as h, p, N, r]m, as well as the operators A and 
Ci defined below, will be often suppressed to simplify notation. 

Problem 3 (Scattering problem, variational form). Given a pair {k,uj), find a function u € Hp^^.{Q) 
such that 

h{u, v) — uP'h{u, v) = p{v), for all v e Hp^^{il). (8) 

By definition, a generalized guided mode is a nonvanishing solution of Problem [3] with p set to 
zero. Such a solution possesses no incident field and therefore satisfies the outgoing Condition [l} If 
(k, uj) is a real pair, then all propagating harmonics in the Fourier expansion Q of the mode must 
vanish and the field is therefore a true guided mode, which falls off exponentially with distance 
from the slab. This can be proved by integration by parts, which yields a balance of incoming and 
outgoing energy flux. In the case of one propagating harmonic, {h,ijo) G P, this means 



\A-,\^ + \B^\^ = \A-^\^ + \B^\\ 



Guided modes with cj 7^ are fundamental in the theory of leaky modes [9l W7\ [TOl [5] and always 
are exponentially growing as 1^1 — t- 00 and decaying in time. The following theorem is proved in 
[T5l Thm. 5.2] and O Thm. 15]. 



Theorem 4 (Generalized modes). Let u he a generalized guided mode with k S M (and Re(u;) > Oj. 
Then Im(cj) < 0; and u — )• as |z| 00 if and only if uj is real. 

It is convenient to write the form a — uP'b as 

h{u, v) — uP'b{u, v) = ci(n, v) + C2{u, v), 

in which ci{u,v) = h{u,v) + b{u,v) and C2{u,v) = — (w^ + l)b{u,v). Both ci and C2 are bounded 
forms in H^^^{Q). If we take V to be a sufficiently small complex neighborhood of the diamond V, 
c\ is coercive for all (k,w) in T)' . These forms are represented by bounded operators C\ and C2 in 
iip\i.(0): 

(Ciu,v)Hpi^^(Q) = ci(n,f). 
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{C2U,v)Hi^^(n) = C2{u,v). 

If we denote by w™^ the unique element of H^^^{^) such that {w^^^ , v) = p{v), the scattering 
problem becomes {Ciu,v) + {C2U,v) = {w^^^,v) for all v G Hp^^^Q), or 

Ciu + C2U = w'''''. (9) 

Because of the coercivity of ci and the compact embedding of Hp^j.{Q) into L^(ri), we have 

Lemma 5. The operator Ci has a bounded inverse and C2 is compact. 

By means of the Fredholm alternative one can demonstrate that, even if a slab admits a guided 
mode for a given real pair {k,u}), the problem of scattering of a plane wave always has a solution. 
Proofs are given in [H Thm. 3.1] and \13\ Thm. 9]; the idea is essentially that plane waves contain 
only propagating harmonics whereas guided modes contain only evanescent harmonics and are 
therefore orthogonal to any plain wave source field. 

Theorem 6. For each {k,uj) G V, the scattering Problem^ with a plane-wave source field has at 
least one solution and the set of solutions is an affine space of finite dimension equal to the dimen- 
sion of the space of generalized guided modes. The far-field behavior of all solutions is identical. 

With the notation 

A{k, uj) = I -\- Ci{k, uj)~^C2{k, w), 
'^ = u and = C7fiu;i'^^ 

equation ^ can be written as 

A{k,u)^{k,(jj) = (I){k,uj), (Scattering problem in operator form) (10) 

in which A is the identity plus a compact operator. A generalized guided mode is a nontrivial 
solution of the homogeneous problem, in which (p is set to zero: 

A{K,u})ip{K,u)) = 0. (Guided mode) (11) 

It can be proved that, if e and /x are large enough in the structure and symmetric in the x 
variable (i.e., about the z-axis normal to the slab), there exists a guided mode at some point (0, ojq) 
in the diamond V [U §5.1]. Such a guided mode is antisymmetric about the z-axis and its existence 
is proven through the decomposition of the operator ^(0, uj) into its action on the subspaces of 
functions that are symmetric or antisymmetric with respect to x. The symmetry of A is broken 
when K is perturbed from zero with the consequent vanishing of the guided mode. 

The frequency uq should be thought of as an embedded eigenvalue of the pseudo-periodic 
Helmholtz operator in the strip {(z,x) : — vr < x < vr} at k = which dissolves into the continuous 
spectrum oj > \K\/y/eoiio as k is perturbed. It is the nonrobust nature of this eigenvalue that is 
responsible for the resonant scattering and transmission anomalies that we study in this paper. 
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2.2 Analyticity 



Analysis of scattering resonance near the parameters {ko,ujq) of a guided mode rests upon the 
analyticity of the operator A. The proof of analyticity is in the Appendix. 

Lemma 7. The operators Ci, C2, and A are analytic with respect to co and k, as long as rj^ = 
u'^eofiQ — (k + m)^ 7^ for all m £ Z, which holds in particular for (k,cj) G V. 

Assume that uj) has a unique and simple eigenvalue i{K, u) contained in a fixed disk centered 
at in the complex A-plane for all (^,0;) in a complex neighborhood of {ko,ujo) G D C M?. (In fact, 
that {ko,ujo) £ is not necessary for the present discussion.) It will be convenient to work with 
i = ci for a nonzero complex constant c to be specified later. 

Given a source field (j){K,uj) that is analytic at {k,q,ujo), consider the scattering problem 



A{k, u})tp 



(12) 



The analyticity of the field is proved in |13| §5.2]. It analytically connects scattering states for 
i{K,(jj) 7^ to generalized guided modes on the dispersion relation £{k,uj) = near (kcwq)- 

Theorem 8. The simple eigenvalue £ is analytic at {kq,ujo), and, for any source field (j) that is 
analytic at [kq^uq), the solution is analytic at {kq,ujq). 

The analytic connection between scattering states and guided modes, introduced in [16], is 
achieved as follows. One analytically resolves the identity operator on Hp^j.(Q) through the Riesz 
projections, 



I = Pl + P2, 

1 
2TTi 



{XI - A{K,uj))-^dX, 



(13) 
(14) 



c 



where C is a sufficiently small positively oriented circle centered at in the complex A-plane. The 
image of Pi is the one-dimensional eigenspace of A{k,lj) corresponding to the eigenvalue £(k,u;) if 
this eigenvalue lies within the circle C. This eigenspace is spanned by the analytic eigenvector 

in which ^(kq, wo) is an eigenvector of A{kq,loq) corresponding to £{kq,loq). The resolution ( 13|[l4 ) 
provides an analytic decomposition of the source field (j) near {kq,ujq) as (j) = aip + (p2, with 



ail) = Pi(j), (t>2 = P2(t). 
Now, letting A2 denote the restriction of A to the range of P2, one observes that the field 

^ = cat/i -|- £A2^(f)2, 

solves Atp = £(f) : 

A^j 



(15) 



(16) 



1 " 




caip 




iatjj 


AP2 











The Riesz projection naturally decomposes the source and solution fields into "resonant" and 
"nonresonant" parts via ( 15|16 ). 
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3 Resonant transmission for a symmetric slab 

From now on, we will assume that the structure is symmetric about the x-axis. Thus, in addition to 
the conditions ([l]), we assume also that e{x, —z) = e{x, z) and —z) = fi{x, z) for all x. We also 
assume that £{k,uj) is a simple (necessarily analytic) eigenvalue in a neighborhood of {kq,u)q) G V 
and that i{Ko,ujo) = 0. 



3.1 The reduced scattering matrix 

For {k, lo) E P, consider the problem of scattering of the field e'^i'-^^+'noz) i^cicient upon the slab on the 
left. By Theorem [gJ a solution exists and the difference of any two solutions is evanescent; in fact, 
near {kq,u)q) the solution is unique if and only ii £{k,uj) / 0. Thus the propagating components of 
the periodic part u of the solution u = e^'^^u are unique, resulting in well-defined complex reflection 
and transmission coefficients R and T, 

u = 6^''°^ + i2e-'^o^ + c" e^^'"^-''-^) for z < -L, 
u = Te*'?"^ + c+e^(™^'+'''"^) for z > L. 

m/O 

Because of the symmetry of the structure with respect to z, an incident field e'^i'^^-voz) fj-Qj^ ^j-^g 
right produces identical reflection and transmission coefficients. Thus the reduced scattering matrix 
for the structure for (k, uj) G V is 



S{k, oj) 



T{k,uj) R{k,uj) 
R{k,uj) T{k,uj) 



which gives the outward propagating components in terms of the inward propagating components 
in the expression Q via S{Aq,B^Y = {A^,BqY. 

In terms of the transmission coefficient T, we define the transmittance to be the fraction of 
energy flux that is transmitted across the slab. The transmittance is the quantitiy |Tp, which lies 
in the interval [0, 1]. 

Let us now take the incident field from the left to be £{K,u})e^^'^^~^^'^^\ which results in a 
refiected field a{K,uj)e''^'^^~^°^^ for z — )• — oo and a transmitted field b{K.,uj)e^^'^^~^^°^^ for z — )• oo, 
with coefficients given by 

a = £R, b = £T. (17) 

By the structural symmetry, an incident field £(K,uj)e'^^'^^~'^°^^ from the right results in a refiected 
field a{K,u})e^^'^^~^^°^^ for z — )■ oo and b{K,uj)e^^'^^^'^'^^^ for z — t- — oo, with coefficients also given 



by (17). The utility of working with a and b is that they are analytic, whereas R and T are not 



analytic at points (k, to) where £ = . 

Lemma 9. The coefficients a{n,uj) and b{K,uj) are analytic in k anduj. 

Proof. The analyticity of the incident field £{K,u})e^^'^^~^''T^^^ implies the analyticity of the source 
field £{k,uj)4>{k,lo) in the equation A{k,uj)xP = £<p and hence, by Theorem [sj also the analyticity 
of the solution field 'ip{K,u) = u{x, z; k,uj) in H^^^{^). The coefficients a(K, w) and b{K,uj) of ip are 
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given by 



^irio{K,uj)L f-2n 

a{K,u}) = / u{x, —L; K,uj)dx, 

27r Jo 

^—irio{K,u})L r2Tr 

b{K,uj) = / u{x, L; K,uj)dx, 

27r Jo 

and since rio{K,oj) is analytic and u i— ?■ j^^ u{x, ±L)dx are bounded linear functionals on Hp^j.{^}), 
both a and b are analytic. □ 

This lemma provides a representation of the scattering matrix as the ratio of analytic functions 
in a complex neighborhood of D, except at points of the dispersion relation i{K, u) = 0, 



S{k, Ijj] 



b{K,uj) a{K,u}) 
a{K,Ljj) b{K,u}) 



(18) 



Assuming that £ = at an isolated point (kq, wq) of 2?, we see that, in a real punctured neighborhood 
of {kojUq), S{k,uj) is a complex-valued real-analytic function. In fact, for real {k,u}), S{K,io) is 
unitary, a standard fact that is shown by integration by parts: 

/W'^l-l' + l"!'- (19) 
\ab + ab = 0. 

This implies, in particular, that three analytic functions vanish at (ko,wo): 

£{ko,ujo) = a{Ko,ujo) = b{KQ,uJo) = 0, 

which is the feature that leads to the sensitive behavior of the transmission and reflection coefficients 
and \a/i\'^ near {ko,ujo)- 
Now, each of a{K,uj) and b{n,uj) is a complex functions of two complex variables, and a 
In this section, we analyze the generic case 

d£ , da , db , , , , , 

9^^°' 

The Weierstrafi Preparation Theorem tells us that the zero-sets of a, b, and i are graphs of analytic 
functions of k near kq. Let k = k — and Co = u) — loq- With the appropriate choice of c in ^ = ci, 
the Theorem yields the following factorizations: 

a{K,uj) = {Cj + rik + r2K^ H )(roe*'^ + r^K + rcjCj + 0{\k\'^ + |a)p)), 

b{K, uj) = {iJj + tiR + t2R^ + ■■■ ){itoe''^ + t^K + + 0(|kP + |tD|^)), (21) 
£(k, uj) = {Cj + Ilk + l2k^ + • • • )(1 + 4k + + 0(|k|2 + IcDp)), 

in which < rg < 1 and either 0<io<lor— l<to<0; the symbols r^, rjj, etc., refer to 
constants. One shows that the same unitary number e^'^ appears in the second factors of both a 
and b by using the second expression in ( 19 ). The zero-set of the first factor of each in each of these 
expressions coincides with the zero set of the corresponding function a, 6, or £ near [ko,loo). 

Under these conditions, one can deduce several properties of the coefficients; see ^16j and [131 
Theorem 10] for proofs. 
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Lemma 10. The following relations hold among the coefficients in the form (21); 

(i) + tl = l, 

(ii) ^1 = n = ti G M, 

(Hi) Im(^2) > 0, 

(iv) ^2 G M r2 = t2 G M ^2 = ^2 = ^2 S M. 

When Ko = 0, the coefficient £1 necessarily vanishes because of symmetry of the dispersion 
relation ^(k, tj) = in k. Whether 7^ can be realized for dielectric slabs remains an open 
problem. Nevertheless, we do not assume £1 = 0. We also assume that Im(^2) > 0, which is 
sufficient to guarantee that {kq,u:q) is an isolated point of i{K,uj) = in D. 

The proof of the following technical lemma is in the Appendix. 

Lemma 11. Under the assumptions (20) and lm(£2) 
satisfied: 

(i) r2 and t2 are distinct real numbers; 

(ii) r2 =t2 M and either £2 = f2 = ^2 or £2 = "^2 = ^2- 



> 0, one of the following alternatives is 



3.2 Resonant transmission 

We now present and prove the main theorem on total reflection and transmission by symmetric 



periodic slabs. The content of the three parts of Theorem 12 can be paraphrased as follows. 



(i) If the coefficients r2 and t2 of the quadratic part of the expansions ( [22p3l ) are distinct 
numbers (case (i) of Lemma 11), then all of the coefficients r„ and t„ of both expansions turn out 
to be real numbers. The consequence of this is that a(K,u}) and b{K,uj) vanish along real-analytic 
curves in D given by 



kq) 



t2{K - 



^0)' - 



(a(K, uj) 
{b{K,uj) -- 



= 0) 
0) 



(22) 
(23) 



Because r2 7^ t2, these curves are distinct and intersect each other tangentially at (kqj'^o) with 
one of them remaining above the other. They give the loci of 100% transmission (|6/^| = 1 when 
a = 0) and 0% transmission (|6/^| = 0). Thus we are presenting a proof of total transmission and 
reflection at two nearby frequencies near loq, whose difference tends to zero as k tends to kq. This 
establishes rigorously the numerically observed transmission spikes in [161 113] . The zero sets of a 
and b are depicted in Fig. [sj The scattering matrix S{k,u!) is not continuous at {kojUJq) because 
|r| = \b/i\ takes on all values between and 1 in each neighborhood of {kq,ujq). However, as a 
function of a; at k = kq, S{hq,u}) is in fact continuous. 



(ii) The result of part (i) of Theorem 12 provides a local representation of the zero sets of a 
and b about (ko,wo) as the graphs of the real analytic functions ( 22|23 ) in V. These locally defined 
functions can be extended to real-analytic functions of k in an interval around kq up until their 
graphs either intersect the boundary of the diamond D or attain an infinite slope. It is not known 
which of these possibilities actually occur in theory. 

(iii) In case (ii) of Lemma 11, we prove that the transmittance |6(k, w)/£(k, a;)p is continuous 
at {hq,u}q). We are aware of no examples of this case but are not able to rule it out theoretically. 
This result tells us that, if there is in fact a resonant transmission spike, then case (i) of the Lemma 
must hold and thus this spike attains the extreme values and 1 by part (i) of Theorem 12 In 
short, either \b{K,uj)/£{K,uj)\'^ is continuous at (ko,<^o) or it attains and 1 in every neighborhood 
of (ko,wo). 
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Theorem 12 (Total transmission and reflection). Given a two-dimensional lossless periodic slab 
satisfying ([T]) that is symmetric about a line parallel to it. Let (kqjUJo) be a wavenumber-frequency 
pair in the region T> of exactly one propagating harmonic at which the slab admits a guided mode, 
that is £(ko,wo) = 0. Assume in addition the generic condition (20) and that Im(£2) > in the 



expansion of £ in (21). Then either the transmittance is continuous at {kqjOJo) or it attains the 



magnitudes of and 1 on two distinct real-analytic curves that intersect quadratically at {kq^ujq). 
Specifically, 

(i) If r2 7^ t2, then rn and t^ are real for all n and both a{KQ,uj) / 1{kq,uj) and 
b{Ko,uj)/£{KQ,uj) can be extended to continuous functions of oj in a real neighborhood of ujq with 
values rge*''' and it^e^^ at ujq. 

(a) Ifr2 7^ t2, let c denote either a orb and letuj = /(k) denote the corresponding function from 
the pair (2^[2§. Then f can be extended to a real analytic function g{K) on an interval (ki,K2) 
containing kq such that the graph of g is in V and for each i = 1,2, the limit g{ni) := lim^^^^. g{K) 
exists and either {Ki,g{Ki)) is on the boundary ofT> or ^{k^, g(Ki)) = 0. 

(Hi) Ifr2 = t2, then \ a/i\ and \ b/i\ can be extended to continuous functions in a real neighborhood 
of {k,q,uo) with values |ro| and |to| ctt {kq,ujq). 



Proof, (i) Assume r2 ^ t2 G M. Recall from Lemma 10 that ri = ti = ii G M. Assuming 
ri, . . . , r„ G M for n > 2, we will show that rn+i G M. For (k, uj) subject to the relation u + iik + 
r2Kp + r^k^ + . . . + r„K" = 0, 



a 
b 



{t2 - r2)K? + (ts - rs)^^ H h (tn - r„)K" + t„+lK 



n+1 



+ • 



roe'^ + 0{\R\) 
itoe'^ + 0{\k\) 



■An-1)_ 



rn+l + rn+2K + 



(t2 - r2) + {ts - r^)R + 



roe'^ + 0(|At|) 
itQC^^ + 0{\k\) 



Because a6 + a6 = (see ([l9|)), it follows that a/6 G iM, and thus 



r„+i + r„+2K + 



{t2 - r2) + (t3 - rz)h + 



roe*^ + 0{k) 
itoe^'y + 0{k) 



G il 



Letting k — )• yields 



Tn+l 



G il 



t2 - r2 ito 

which implies that r„+i G M. We conclude by induction that r„ G M for all n > 1. The proof that 
t„ G M for all n is analogous. 

To prove the second statement, one sets k = n — kq = and observes that the ratios a/i and 
b/i have limiting values of rge*'^ and it^e^'^ , respectively, as — t- 0, or a; — t- ojq. 

(ii) Define the set 

G := {g : (k_, — )-R | kq G (k-, g is real-analytic, a{K, g{K)) = 0, T{g) G V}, 

in which T{g) is the graph of g, and define the numbers 

Ki := inf{K_ \ g gQ}, 
K2 := sup{k4. \ g eG}. 



By virtue of the function (22), which belongs to Q, ki < kq < K2. Standard arguments show 



that any two functions from Q coincide on the intersection of their domains, and one obtains 
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thereby a maximal extension g £ Q oi (22) with domain (ki,K2). We now show that \iTS\i^y^i^^ di^) 
exists. Set uj- = liminff^y^i^^ g{K) and ijJ+ = hmsup^^^^ Because of the continuity of g, 

the segment (^2, ['^-7'^+]) in ^ is in the closure of the graph of g, on which a vanishes. Thus 
0(^2,0;) = OVw G L1J4.]. Moreover, for each uj E there is a sequence of points 

with Uj K2 and a{Kj,uj) = from which we infer that d"a/dK,'^{K2,uj) = OVn G N and hence 
that d^~^^a/doj'^dK.'^{K2,uj) = \/m,n G N and G (a;_, If {uj^,uj^) is nonempty, then 
a must vanish in D, which is untenable in view of the assumption that da/du}{Ko,uJo) 7^ 0. This 
proves that =uj^ so that 0)2 '■= limK/'K2 di^^) exists. If |«;|/^eo/^o < W2 < (1 — l'^l)/\A^oW) a-^id 
da/du}{K2,uJ2) 7^ 0, the implicit function theorem provides an element of Q with k+ > K2, which is 
not compatible with the definition of K2. Analogous arguments apply to the endpoint ki and to 
the function b. 

(iii) If r2 = t2, then by Lemma [111 t2 = £2 or t2 = £2- Keeping in mind that Im(£2) > and 
£1 G M and restricting to (k, u) G 



1)2 



lim 

K.aj— >0 



W + ilk + t2K^ + tsK^ + 



lim 

K,aj->-o a; + liR, + ^2^^ + ^3^^ + 



itoe'^ + 0{\Co\ + \k\) 



: lim |to| 



a) + ^IK + t2K^ 



UJ + ^IK + . 



Whether t2 is equal to £2 or £2, the second factor of the last expression is equal to unity, and we 
obtain limK,(i,-s>o \b/£\ = \to\- Similarly, one shows that limK_^_j.o = |?"o|- D 



3.3 Discussion of assumptions and conditions 

A number of assumptions and conditions concerning the slab structure, the nature of the dispersion 
relation i{K,uj) = 0, and the coefficients in the Weierstrafi representations of i{K,uj), a{K,uj), and 
b{K,uj) have played roles in our analysis. We take a moment to discuss their relevance and the 
regimes of their validity. 

Existence of embedded guided-mode frequencies. Real (^,0;) values at which a periodic slab 
supports a guided mode are typically located below the light cone w^eo/^o = for the ambient 
medium in the first Brillouin zone, that is below the diamond D. The reason is that the light 
cone marks the bottom of the continuous spectrum for the K-periodic Helmholtz operator in the 
strip {{x,z) : \x\ < vr} jl] and thus any guided mode frequency above the light cone corresponds 
to an embedded eigenvalue and is therefore unstable with respect to generic perturbations of the 
operator. Embedded guided-mode frequencies arise under special conditions; in particular, struc- 
tures symmetric in x for which e and ^ are sufficiently large admit a positive finite number of such 
frequencies in 2? at k = [Ij, which are antisymmetric in x. 

Symmetry. The assumption of symmetry of the structure about a line parallel to it is crucial 
to the proof of total reflection and transmission. The condition is used in part (i) of the proof of 



Theorem 12 when the fact that a/b is imaginary is invoked. This is a technical step in the proof 



and does not provide much of an enlightening explanation of the role of symmetry. At least one 
can say that the functions ( 22|23 ), which are in general complex- valued, become real- valued when 
the structure is subject to a certain condition of symmetry. 

Simplicity of the eigenvalue £. The frequencies of these antisymmetric modes are eigenvalues of 
a modified scattering problem in which the one imaginary multiplier — ir/o = — ^|%| in the definition 
of the Dirichlet-to-Neumann operator A'^ ([g]) is replaced by zero. This results in a symmetric 
modification hr{u,v) of the non-symmetric form h{u,v) defined in ([T]). For k = 0, the associated 
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operator is invariant on the subspace of H^^^{^) consisting of antisymmetric functions in x. Thus it 
can be viewed as a Dirichlet operator in the strip {{x, z) : < x < vr}, and the smallest eigenvalue 
is therefore simple. 



The generic condition (20) prevails in numerical computations. The nonvanishing of di/du 
actually implies the nonvanishing of either da/duj or db/doj. If one of the latter vanishes, the 
limiting value of the transmitted energy at {kqjUq) is either 1 or 0; we will consider this case in the 
next section. Verification of any of these conditions for a particular structure is done numerically 
(or by experiment). 

Nonrobustness. Guided modes at embedded frequencies tend to be nonrobust under perturba- 
tions of K, which means that the point (kq, wq) is isolated in the domain D. Typically, one computes 



that Im(^2) is strictly positive in (21), which guarantees that (ko,wo) is an isolated point in the 
intersection of D with the dispersion curve £{k,uj) = 0. So far, we do not have an analytic proof of 
nonrobustness for any particular structure. In fact, one can show that real dispersion relations for 
robust guided modes do exist within the continuous spectrum if the structure possesses a smaller 
period than that of the pseudo-period of the guided modes. 

Traveling modes and ii 0. A guided mode at k = is in fact a standing wave, and because of 
the symmetry of £{k,u}) in k, necessarily ii = 0. In this case, the anomaly is not linearly detuned 
as K is perturbed from but simply widens at most quadratically in k. As far as we know, guided 
modes that are not robust under a perturbation of k, have not been demonstrated for kq ^ or 
for structures that are not symmetric in x, although we believe that they should exist in slabs 
with anisotropic components. Guided modes with kq ^ have been proved in a discrete system 
in which the ambient space is a two-dimensional lattice and the waveguide is a one-dimensional 
periodic lattice coupled along a line to the 2D one and for asymmetric scatterers in another 
discrete model [14J. In fact, both of these systems are amenable to the analysis in this section and 
are subject to Theorem 12 with £i not necessarily zero. Compare Figures |4] and [5] 



The alternative r2 = t2 ^ M. oi Lemma 11 Under this condition. Theorem |12| states that the 
transmittance exhibits no spike; in fact, b/£ is continuous at {kq,ujo). Although we are not able to 
rule out this alternative, we are aware of no examples in which it does occur numerically. In any 



case. Theorem 12 guarantees that, if the transmittance is discontinuous at (kq,ojo), it must attain 



the values of and 1 along the real-analytic curves ( [22p3l ) 



4 Nongeneric resonant transmission 

The anomalies we have analyzed so far possess "background" reflection and transmission values of 
|ro| and |io|) both of which lie strictly between and 1. These are the limiting values of the ratios 
of reflected or transmitted energy to the incident energy at k = kq and a; — )• cjq (see part (i) of 
Theorem 12). When one of these values reaches or 1 (we always have |rop + |toP = 1)) the form 



for a or 6 in (21 ) must be modified so that the order of u in its first factor exceeds one. This results 
in a transmission anomaly possessing a single sharp dip or peak near the frequency of the guided 
mode — the peak-dip feature reduces to an upright or inverted Lorentzian shape (Fig. [6|[7|. 

We also discuss what happens when i is degenerate at (kq,ujq). We still assume that ^ is a 
simple eigenvalue of the operator A{k,uj) but allow the first partial derivative of i with respect to 
u to vanish at {kq,ujo). Assuming that the second partial is nonzero, two double-spiked anomalies 
instead of one emanate from the guided mode frequency upon perturbation of k (Fig. [sj [o]). We 
do not have a rigorous demonstration of a case of degeneracy, but we have observed numerically a 
pair of anomalies emerging from a single guided-mode frequency for a slab consisting of two infinite 
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rows of identical dielectric rods side by side. 

Wc shall omit the proofs, which are technical. Resonant transmission in nongencric cases is 
worthy of an in-depth investigation in its own right, and our aim here is to provide a view of some 
of the phenomena that can theoretically occur. 



4.1 Total background reflection or transmission 

Let us assume that dl/duj 7^ at (kq, then only one of the functions a and h can be degenerate 
at {kq,ujq). 

Proposition 13. Suppose that at {kq^ujq) G T>, 

di 

^(/«o,wo) = 0, and — (^0,^0)7^0. 

Then 

da \ / r, / 

— {ko,loo)^0 or — (ko,c^o) 7^ 0. 

Proof. By conservation of energy |£p = |ap + |6p for {k,u) G V, we have a(/io,wo) = and 
6(ko, Wo) = and therefore the representations i{Ko,Lo) = {u — uJo)gi{uj), a{0,uj) = {cj — iuo)'^g2{co), 
and 6(0,0;) = {lo — a;o)"53(oj), with gi analytic and nonzero at ujq. This is consistent with \£\'^ = 
|ap + |6p only if m = 1, that is da/duj 7^ 0, or n = 1, that is dh/doj 7^ 0, at (ko,u;o). □ 

Let us consider the case of 100% background transmission, that is dl/duj, d'^ajduj^ ^ and dh/du 
are nonzero and that da/doj = at {kq^uq). The Weierstrafi Preparation Theorem yields the 
following factorizations: 



w) = (cD + + i2h'^ + • • • ) (1 + 0{\k\ + lo-D), 
a{K,uj) = {Co^ +Oja^{R) + a^{R)) {roe'" + 0{\R\ + \Cj\)), 
b{K, u) = {oj + tiK + t2K^ + ■■■) {toe'^ + 0{\R\ + |w|)), 



where ro, to > 0. 

Let us also suppose that Cj"^ + Cja^{R) + a^{R) has distinct roots at k = 0, so that it can be 
factored analytically, 

a={u + r^^h + + + rfh + rf^R^ + • • • ) (roe'" + 0{\R\ + |a;|)) . 

Assuming that the coefficients are real, maximal transmission is achieved now at two points 
that move apart linearly m. R = k — kq. 

Lemma 14. The coefficients of the expansions satisfy the following properties. 

(i) to = l,ti=£ie R, Im(^2) > 0; 

(ii) e*^ = ±ie*"; 

(Hi) r^^ + r^^\ r^^V^^'', and (r^^'' —ii){r^i'^ <3ire real-valued. Therefore, either r\^\ r^^ are both 
real or they are conjugate complex numbers. 

One of the following alternatives holds. 



(i) If (rj^^ - h){rf^ - £i) + 0, then t^ = Re{i2) and |Im(£2)P = 

(ii) If {r\^^ - £i)(rf ^ - 4) = 0, then Re(£2) = Re(t2) and |Im(£2)| = |Im(t2 
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We obtain the factorizations 



a 



+ £ik + e2K^ + ---){1 + 0{\k\ + \u\)), 
UJ + r^^^R + ^^^^2 + ...^[cj + rf^K + rf^k^ + • • • ) (roe'^ + 0{\h\ + |cD|)) , (24) 
6 = (tD + ^iK + ^2^2 + • • • ) (±ie'^ + 0(|k| + . 

If all the coefficients t„ are real- valued, which is guaranteed by the first alternative (see Theorem 
15), then the transmission vanishes along the real-analytic curve 

UJ = UJQ — £i{k — Kq) — t2{K — Kq)^ — ■ ■ ■ . (5 = 0) (25) 
If and r^^^ are real-valued then the transmission achieves 100% along two real-analytic curves 



UJ 



ujQ-rX'in- Ko)-r)2'{K- kqY i = 1,2. (a = 0) (26) 



The two frequencies of total transmission and the one frequency of total reflection come together 
at Wo as K tends to kq. The frequencies of maximal transmission are difficult to detect in Fig. [6] 
and [7] because the transmission function is nearly constant near these points. A magnified view is 
shown on the right of Fig. [7| The proof of the following theorem is similar to that of Theorem 12 
We believe that ti must lie between r^^'' and r[^^ but do not have a proof of this. 

Theorem 15. Suppose that e and fj, are symmetric in z and that i(KQ,ujQ) = at (kq,u}q) £ T>. Let 

sb' ^'"''^ §u} nonzero and ^ = at {kq,ujo). Suppose in addition that (r^ —ii) / 0. 

Then^ 

(i) tn is real for all n and therefore the coefficient b of the transmitted field vanishes on the 



real-analytic curve in a neighborhood of {kq,ujq) given by (25); 



(a) if r^i^ and r^^ are distinct real numbers, then rn^ and rn^ are real for all n and therefore 



the coefficient a of the reflected field vanishes on the real analytic curves given by (26). 
4.2 Multiple anomalies 

If both I and dl/duj vanish at {kq,ujq) G X>, then again by virtue of the equation |£p = |op + |6p 
for real (k,(x>), both a and b and their (^-derivatives vanish at {hiq,U}q). We treat the case in which 

dl da db 

^ = 0, — =0, — = 0; 

OUJ OUJ OUJ (cyrj\ 

The zero loci are given locally by the roots of a quadratic function in uj with coefficients that are 
analytic in k and vanish at k = 0, 

£{k,uj) = {u^ + LoX^{k) + A°(k)) (1 + 0{\k\ + \Cd\)), 
a{K,uj) = {Cj^ + Cba^ik) + a°{k)) (roe*^ + 0{\k\ + |a;|)), 
b{K,uj) = {u^+cbl3^{k) + I3\k)) (it^e'^ + 0{\k\ + \uj\)), 

in which A*(0) = a*(0) = /3^(0) = 0, < ro < 1, and to is real with < \to\ < 1. The ratio of 
the leading coefficients in the factors of a and b is imaginary because of the special form of the 
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scattering matrix (18) coming from the symmetry of the structure in z. Let us assume again that 



the left-hand factors have distinct roots so that they can be factored analyticahy. Here again, we 
are able to show that just two real coefficients and for the linear terms k suffice, as long as 
they are distinct, 

t{K,uj) = (^w + 4^)k + 4^)k2^---) (^ti + 4^)^ + 4^^^^ + •••) {'^ + 0{\R\ + \Cj\)), 

a{K, uj)=(uj + 4^^^ + 4^^/^^ + ■■■') {cj + ifk + rf-K^ + • • • ) {r^i^ + r«K + r^tl; + • • • ) , (28) 

6(/t, uj) = (cj^ i^l^R + 4^^^^ + . . . ^ (^cj + if>k + 4^^^2 + • • • ^ (itoe*'^ + t^^K + t<i6D + • • • ) . 

Lemma 16. The coefficients 1^^^ and i^f^ are real and the coefficients 4^^ md i^2^ have nonnegative 
imaginary parts. Moreover, i/ Im(^2^)) > and Im(£2^)) > 0, then, for each i G {1,2}, either 

(i) (i) 

r2 and distinct real numbers or they are equal and not real. 

Subject to structural symmetry with respect to z, one can prove a theorem analogous to those 
for the previous cases. It says in particular that, if the first alternative in the last statement of the 
lemma holds for i = 1 and i = 2, then two peak-dip anomalies emerge from ujq as k is perturbed 
from kq, each of which attains the values and 1 along real-analytic curves passing through (kq, ujq). 

Theorem 17. Suppose that Im(4*'') > for all i e {1,2}. 

(i) If r2^ 7^ t2'' with i G {1, 2}, then rn^ and t^ are real for all n. 

(a) If r2^ = for all i G {1,2}, then aji and h/l are continuous at {k,q,uio). 

5 Transmission graphs 

We demonstrate resonant transmission anomalies with various choices of the coefficients in the 
expansions ( |21p4p8| ) of I, a, and h. We graph the transmittance 



\T{k,u) 



|a|2 + |6|2^ 



(29) 



retaining terms up to quadratic order m. h = k — kq m. the factors that vanish at {kq,ujq) and only 
the constant in the nonzero factor. An error estimate of 0(1^1 + uP') for these approximations in 
the generic case is proved in fll'j Thm. 16]. First we demonstrate the nondegenerate case of Section 
3.2, in which r2 and t2 are distinct real numbers (Fig. |4| [5]). For k = kq {k = {)) the anomaly is 
absent; it widens quadratically in k, namely as \t2—r2\i^. If H-i 7^ 0, as in Fig. [sj then the anomaly 
is detuned from ujq at a linear rate in R, whereas it widens only quadratically. 

Figures [6] and [7] show the degenerate case in which the anomaly is a single dip descending to 
zero from a background of full transmission or a single peak rising to unity from a background of no 



transmission (Section 4.1 ). In the former case, for example, full transmission is actually achieved at 
precisely two frequencies near uj = ojq (o) = 0); this can be seen clearly in the magnified, right-hand 
image of Fig. [7j 

When the first derivative with respect to a; of all three functions i, a, and h vanishes at (kcwq), 
an anomaly characterized by a pair of peak-dip features can emanate from wq, as discussed in 



Section [42| Figures|8]and|9]illustrate two possible choices of constants in the Weierstrafi expansions. 
At this point there are no rigorous results on which sets of constants are admissible or prohibited 
by periodic slabs or any other open waveguide system. 

The vertical line in all graphs shows the location of a) = 0, or cj = wq- 
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Figure 4: jrp as a function of Cj for k = 0, ±0 01 ±0 .02, ±0.03. The generic conditions ^, f^, |^ 7^ 
are satisfied at the bound-state pair [kq^loq). In ( 22|23 ), £1 = so that there is no hnear detuning of the 
anomaly with k. Left: < ^2 = 1 < ''2 = 2 so that the peak is to the left of the dip and both are to the left 
of ujQ. Right: r2 ~ —2 < < i2 = 1- In both graphs, tq = 0.6, tf, = 0.8. The transmission is symmetric 
in K, and the curve without an anomaly (the horizontal curve) is the transmission graph for k = 0. 
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Figure 5: |T|2 as a function of w for k = 0, ±0.00 3, ±0.0 06, ±0.009. The generic conditions 
are satisfied at the bound-state pair (kqjWo). In ( 22|23 1, l\ — 0.9 7^ 0, so the anomaly is detuned from 
a; = wq (w = 0) in a linear manner in k. The coefficients r2 = 2 and ^2 = 1 of k? are distinct real numbers, 
and (ro, t^) = (0.6, 0.8). The graph for k = is horizontal at jTp k. 0.64. 
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Figure 6: 
occurs when 

and < < rx 



|rp as a function of uj for k = 0, ±0.01, ±0.02, ±0.03. Left: Full background transmission 
£ ^ 0' Is = 0' and £ ^ at («o,^o). In g, < rf^ = 0.2 < ii = £1 = 2 < rf^ = 4, 

£^0, IS^O, and£=Oat(^o,c.o) 



(7, 7, 0.1), and = 0.6. For k 0, \T\^ « 1. Right: 



Ix < t^^\ For k = 0, |r|2 « 0. 
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Figure 7: |Tp as a function of u for k = 0, ±0.01, ±0.02, ±0.03. Left: Full background transmission occurs 



when 



(1) .(2) 



^ 0, - 0, and £ ^ at (ko,c^o). In (24), 



.(1) - 



= -0.04 < h = £i = < = 0.06; 

t^) = (0.6,-1,1,1). For K = 0, |r|2 « 1. Right: Magnification of the graphs for k = ±0.01, 



bringing into view the frequencies of total transmission. 




-0.0030 -0.0028 -0.0026 -0.0024 -0.0022 -0.0020 -O.OOIS 



Figure 8: Left: jTp as a function of w for k = 0, ±0.003, ±0.006, ±0.009. The partial derivatives of i, a, 



and b all vanish at (kq, 1^0)7 whereas their second derivatives are nonzero. In (28l, (^i^'',^!^'') = (0.7,0. 



(1) 



(ro,to) - (0.6,0.8), r^'^ = 2<t 
above |T|2 = 0.6. Right: k = 0.003. 



8 and t 



(2) 



4 < r. 



(2) _ 



5. For k = 0, the graph is horizontal, just 





-0.0040 -0.0035 -0.0030 -0.0025 -0.0020 -0.0015 -0.0010 



Figure 9: Left: |r|2 as a function of w for k = 0, ±0.003, ±0.006, ±0.009. The partial derivatives of i, a, 
and b all vanish at {kq,ujq), whereas their second derivatives are nonzero. In (28l, {£^i\£'"i^) — (0.7,0.8), 
(0.6,0.8), r^^^ = 2 < = 8 and = 4 < t 



(2) 



6. For K = 0, the graph is horizontal, just 



above \T\^ = 0.6. Right: k = 0.003. 
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A Appendix 

Proof of Lemma^ To prove Ci is analytic with respect to oj at (kq, ojo)-, we let u = ujq + Aw and 
show that 

Ci(ko, Wo + Aw) - Ci(Ko,a;o) 



lim 

Aw-5>0 



Aw 







in operator norm, where Ci(^(ko,wo) : H^^^{^1) — )• //pgj,(r2) is defined by 



-E 



-iep^owo 



The partial derivatives of r/m with respect to w are 



and the joint analyticity of ry.m at (kq, wq) implies that 

T]m{K0,U)0 + Aw) = r/m(Ko, Wo) + ?7mL^(Ko, wo)Aw + i?(^^(Aw) 



?/m(Ko, Wo) + 



with the remainder term 



fAw)2 



27ri 



eo^o^o 
V^eo^owg - (m + Ko)2 

^m(s) 

Co (S - w)(s - Wo)2 



Aw + i?(^)(Aw) 



in which Cq is the circle in the complex plane centered at wq with radius rg and 2|Aw| < tq. For 
any s £ Co, |s — w| > |ro — |Aw|| > tq — |Aw| > ro/2. We estimate that for all m ^ 0, 



i?(^)(Aw) 



< 



|Aw|' 
2tt 



2vrro 



suPseCo{l^m(s)|} ^ sup,gCo{l^m(s)|}(|Aw|)2 
(ro - |Aw|)rf 



ro/2 • ro 



In the last expression, if m = 0, |?/m(s)| < C for some C > 0. If m 7^ and eo/UoWp — (m + k)^ < 0, 
|??m('S)| = ^y {m + kY — eo/^ow^ < (m+ft) < 2m < Cm for some C; if m 7^ and eo/^oWQ — ("1-+^)^ > 
0, we have |r?m(s)| < \/ eofio^o < C < C*"!. To summarize, Ir/mCs)] < C(m + 1). 
For some constants C'C" > 0, 



Ci(ko,wo + Aw) - Ci(«;o, Wo) 



C 



luj 



n, V 



m 



(r?m(K;o, Wo + Aw) - r/„(Ko,wo)) 



eo/^owo 



ii^^^(Aw) 



Aw 



Aw 

< — Aw 



V^eo^oWo - (m. + ko)^ 

+ l)UmVm 



UmVm 



< C'l Awl ||li|| TTl 



//i/2(r)FII/fi/2(r) 



< C"| Awl llnll iji 
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and 



lim 



C 



= lim sup 



Alo 

Ci{ko,U!0 + Au!)—Ci{ko,U!o) _ ^ 



\u 



= 0. 



and hence Ci is analytic with respect to uj. 

Now we show that the operator Ci is analytic in k. We prove the limit 

Ci(ko + AkjWo) - Ci{ko,u;o) 



lim 

Ak-5>0 



in operator norm, where Ci^ is defined by 

1 . 



Ak 



-C 



Ik 



= 



(CikU, v) = I - [lUVx - lUx 



The partial derivatives of r/m with respect to k are 



-{m + k) 



and 



??m('«o + A«;,wo)-'?m(«;o,'^o) = (m + kq) AAt + r4^)(AAt), 



with |ri^'*(AK)| < D{m + 1)\Ak\'^ for some constant D > and |A«;| sufficiently small. So 
((Ci(ko + AK,a;o)'u,w) - {Ci{kq,u;o)u,v)) 

= / -[V + ^(Ko + Ak)]u ■ [V - i{Ko + Ak)]v - / -(V + iKo)^* • (V - zko)^^ 

+ ^ [VmiKQ + Ak, Wo) — IJmil^O, ^o)] 



fJ-O 



j —[AK{iuVx — iUxV + 2k,quv) + (Ak)^uu] + 



{m + ko) 



and 



/ {Ci{kq + Ak,u:q)u,v) - {Ci{kq,uq)u,v)\ 

\ a;^ ) - 



y'eo/Uoo;^ - {m + kqY 
v) 



A/. + r(?) 



f —UvAk + y V^UmVm < |A/c| 

< D'\Ak\ (||w||Hpi„(n)lbllHpi,,(Q) + ll«ll/fi/2(r)ll 
for some constant D',D" > 0. Therefore 

C\{ko + A/i, Wo) - Ci{kq,ujo) 



-uv 



+ 



Mo|Ak| 

<L>"|AK|||n||^i^^(f^)||^;||^^i^^(n) 



lim 



Ak 



= lim sup 



Ik 



(Cl (KO+AK;,U)o)M,tl) — (Cl (KO,U)o)M,'i;) 

Ak 



= 



21 



in the operator norm and so Ci is analytic with respect to k. 

To prove the analyticity of C2 with respect to to at (ko,wo)i we define an operator C2(^(kO)'^o) 

by 

{C2uU,v)h^^^(P_) = -2ujQ I euv 



and we have 



C2(Ko,a;o + Ao;) - 02(^0,^0) 
Aw 

ojI + 2a;oAa; + Aw^ - 



2aj 



Aw 



+ 2wo 



-Aw 



euv. 



As Aw — )• 0, this tends to 0, and thus C2 is analytic with respect to w. 

The operator C2 does not depend upon n. Because Ci is an analytic automorphism, it has an 
analytic inverse and, hence, A is analytic. □ 

Proof of Lemma\ll\ We first compare the coefficients in |£p — |ap — |6p = using the expansions 
(21) and keeping in mind that ^1 = ri = ti. The coefficients of k?, k?Cj, koj"^ and w^ are 

h [(^2 + I2) + 4(4 + 4) - rl{r2 + f2) - 4ro(r^e-^^ + f^.e^T) - tl{t2 + ta) + i^tolt^e"*^ - t^e*^)] , 



(4 + I2) + 2iiih + 4) + ilii^ + 4) - rg(r2 + f2) - 24ro(r^,e-*^ + r^e'^) 
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(4 + 4) + 24(4 + 4) - ro(r^.e-*^ + f;,e^^)-24ro(r<^e-^^ + r^e"^') 

+ itoitf^e-'^ - i-J^) + 2i4io(4e-*^ - 4e'^), 



and 



(4 + 4) - ro(r<ie-^^ + f^e^^) + ito(4e-'^ - 4e'^). 



Define the real quantities A = r^e + r^e*''', S = z(tKe — tae^^), C = r^^e ^'^ + rjjC*''', D 



«(4e-*^ - 4e*'^). Since Re(|£|2 
and w^ are all 0, that is. 



0, the real parts of the coefficients of k^, k^uj, RuP' 



£i[2Re(/2) + 2^iRe(4) - 2rgRe(r2) - IiTqA - 2tgRe(t2) + htoB] = 0, 
2Re(4) + 2^iRe(/s) + 2^2j^e(4) - 2r2Re(r2) - 2hroA+ 

- ijroC - 2tgRe(t2) + 2iitoB + iltoD = 0, 
2Re(4) + 4£iRe(4) - tqA - 2iiroC + toB + 24to-D = 0, 
2Re(4) - roC + toD = 0. 



This linear system can be reduced to 



/O 0\ 



ro to 

\0 ro to J 



B 
C 



/4[2Re(4) - 2r2Re(r2) - 2t2Re(f2)]\ 
2Re(4) - 2r2Re(r2) - 2t2Re(t2) 
2Re(4) 

V 2Re(4) / 



which implies Re(4) — r-QRe(r2) — tQRe(t2) = 0. 
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The coefficients of /c^, k'^Co, kCj'^, and in ah are 

-iroto{r2 + h) + 1h [r^t-J^ - it^v-^e-'^) + 1\ [r^t^e'^ - itor^e-'"') , 
{rohe'^ - itor-^e-'-') + 2h {rotc^e'^ - itova^e-'^) , 

Since ab is purely imaginary, so are these coefficients, and from the second and third, we obtain 

Along the curve {(k, w) € : a) + £ik = 0}, the coefficients of in — |ap — |6p and in ab 
are |£2p - ?'ol''2p - iol*2p and -iroritoh- This yields |£2p = ^^ol^2p + iol*2p and r2t2 G M. 

The relations r2 + t2 G 1^ and r2t2 G imply that either r2,t2 G M or r2 = t2- Because, by 



Lemma 10 ^2 G ^^2 = ^2 G M and because we assume Im(^2) > 0, the numbers r2,42 cannot 
be identical real numbers. In the second case with r2 = t2 from Re(^2) = rgRe(r2) + toRe(t2), 



2 



K2I = ^okal + *ol*2| , and + = 1, we find that Re(^2) = Re(r2) = Re(t2) and |Im(^ 
|Im(r2)| = |Im(t2)|. □ 
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